The response of mechanical and electronic properties of graphane to the elastic strain 
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Based on first-principles calculations, we resent a method to reveal the elastic properties of recently 
synthesized monolayer hydrocarbon, graphane. The in-plane stiffness and Poisson's ratio values are 
found to be smaller than those of graphene, and its yielding strain decreases in the presence of 
various vacancy defects and also at high ambient temperature. We also found that the band gap 
can be strongly modified by applied strain in the elastic range. 



Two dimensional (2D) monolayer honeycomb struc- 
tures of graphenepJ2l BNP and silicorP offer remark- 
able properties and are promising materials for future 
applications. Honeycomb structure of graphene with 
sp 2 bonding underlies the unusual mechanical proper- 
ties providing very high in-plane strength. Graphene 
and its rolled up forms, carbon nanotubes are among the 
strongest and stiffest materials yet discovered in terms 
of tensile strength and elastic modulusP^ Graphane, 
another member of honeycomb structures was theoret- 
ically predicted-^ and recently synthesized by exposing 
graphene to hydrogen plasma discharge P Here each car- 
bon atom being bonded to one hydrogen atom is pulled 
out from the graphene plane and hence whole structure 
is buckled. Instead of being a semimetal like graphene, 
graphane is a wide band gap semiconductor and can 
attain permanent magnetic moment through hydrogen 
vacancies P 

In this work, we revealed the relevant elastic constants 
of graphane using strain energy calculations in the har- 
monic elastic deformation range and compared them with 
those calculated for other honeycomb structures. We also 
found that in the presence of hydrogen vacancy and car- 
bon+hydrogen divacancy, its yielding occurs at smaller 
strains. Furthermore, its band gap first increases then 
decreases steadily with the increasing applied strain. We 
believe that our predictions are relevant for the current 
research focused on the elec tronic properties of honey- 
comb structures under strainP^E] 

First-principles plane wave calculations are carried out 
within density functional theory using PAW potentials.^ 
The exchange correlation potential is approximated by 
generalized gradient approximation (GGA) using PW91 
functional. A plane-wave basis set with kinetic energy 
cutoff of 450 eV is used. All atomic positions and lattice 
constants are optimized by using the conjugate gradient 
method, where the total energy and atomic forces are 
minimized. Interactions between adjacent graphane lay- 
ers in supercell geometry is hindered by a large spacing of 
~10 A. To correct the energy bands and band gap values 
obtained by GGA, frequency-dependent GqWq calcula- 
tions are carried out. Go Wo corrections are obtained by 
using (12x12x1) k-points in the Brillouin zone (BZ), 400 
eV cut-off potential for Go Wo, 160 bands and 64 fre- 
quency grid points. All numerical calculations are per- 




Figure 1: (Color Online) (a) Schematic representation of the 
atomic structure of graphene and graphane. (b) (8x4) rect- 
angular supercell configuration of the system containing 128 
C-H pairs used for the calculation of the elastic constants. 
a x and a y are the lattice constants of the supercell in x- and 
^-directions. Shaded region is the smallest unit cell, (c) The 
mesh of data points (a x , a y ) used for the total energy calcula- 
tions. The units are given in Angstroms, (d) The 3D plot of 
a x ,a y and corresponding total energy values. The red balls 
are actual points and the lines are the fitted formula. 



formed by using VASP package! 13 * 14 * 

The graphene has a 2D hexagonal unit cell with a lat- 
tice constant of a = 2.47 A. The C-C bond length is 
<i=1.42 A and all atoms lie in the same plane. Upon hy- 
drogenation, the lattice constant increases to 2.54 A and 
d increases to 1.53 A. Moreover, C-H bonds are 1.11 A 
and the amount of buckling between the alternating car- 
bon atoms in a hexagon is 0.46 A. Atomic configuration 
of graphene and graphane structures are shown in Fig. [T] 
(a). 

The elastic properties of homogeneous and isotropic 



2 




0.00 0.10 0.20 0J0 0.40 



2> 
2 20 


■(c) /^y 




60 10 
8 
















40 6 








20 ' 


t. s 
iJ 5 






2 


II 











O.ttO 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.1 

Strain s [Ac/c ] 

Figure 2: (Color Online) Two dimensional graphane under 
uniform expansion. (a) Initial atomic configuration in a 
(10x10) supercell treated with periodic boundary condition, 
(b) The variation of strain energy Es and its derivative. The 
orange/shaded region indicating the plastic range. Strains 
corresponding to two critical points in the elastic range are 
labeled as e ci and e C2 . (c) Similar to (b) for a single H-vacancy 
in a (10x10) supercell. (d) For C+H-divacancy in a (10x10) 
supercell. 



materials can be represented by two independent con- 
stants, Young's modulus Y and Poisson's ratio v. Since 
the thickness of a monolayer structure h is ambiguous, 
the in-plane stiffness C is a better measure of the strength 
rather than Young's modulus. Defining Aq as the equi- 
librium area of the system, the in-plane stiffness can be 
given as, C — ~^{ d d % s ), where Es is the strain energy 
calculated by subtracting the total energy of the strained 
system from the equilibrium total energy and e is the uni- 
axial strain (e = Act/a, a being the lattice constant). The 
Poisson's ratio which is the ratio of the transverse strain 
to the axial strain can be defined straightforwardly as 

^ ~£trans / taxial • 

For calculation of elastic constants of graphane, we 



consider large supercell comprising 32 rectangular unit 
cells (8x4). The calculations are also repeated in (2x1), 
(4x2) and (6x3) supercells and the obtained results are 
almost identical, since no reconstructions are observed in 
the system. Fig. [I] (b) shows the supercell used in the 
calculations. a x and a y are the lattice constants of the 
supercell in x- and y-directions in any strain condition. 
In the harmonic region, a's are varied with the strain 
values between +/- 0.02. A grid data (a x ;a y ) containing 
225 points is obtained as shown in Fig. [I] (c). For each 
grid point, the corresponding supercell is fully optimized 
and its total energy is calculated as shown in Fig. [I] (d) . 
By using the least squares method, the data is fitted to 
the formula; Es — a\e x + a2e 2 + a^,t x e y ; where e x and 
t y are the small strains along x- and y-directions in the 
harmonic region. As a result of isotropy in the honey- 
comb symmetry, a± is equal to 02. The same equation 
can be obtained from elasticity matrixP^ in terms of elas- 
tic stiffness constants, namely di = 1X2 = (h- Aq/2) ■ G\\ 
; 0,3 = (h ■ Aq) ■ Ci2- Hence one obtains Poisson's ratio 
v which is equal to C12/C11 = a3/2ai. Similarly, the 
in-plane stiffness, C = h ■ C\\ ■ (1 — (C11/C12) 2 ) = (2a x - 
(a3) 2 /2ai)/(Ao). The calculated values of C by using the 
present method for graphane, graphene, BN, Si and SiC 
2D honeycomb structures are, respectively, 243, 335, 267, 
62 and 166 (J/m 2 ). Also the calculated Poisson's ratios 
are 0.07, 0.16, 0.21, 0.30 and 0.29. Our calculated value 
of the in-plane stiffness of graphene is in good agreement 
with the experimental value^of 340±50 (N/m) and jus- 
tifies the reliability of our method. As seen from the cal- 
culated values, the change of the bonding type from sp 2 
to sp 3 and buckling of the atoms in graphane structure 
makes it 27% less stiffer than graphene. This difference 
can be used to distinguish graphene and graphane ma- 
terials. Also the Poisson's ratio of graphane is almost 
half of the Poisson's ratio of graphene, since the buckled 
structure of graphane reduces the transverse contraction. 
Note that depending on their types and concentrations 
the defects can alter the above elastic constants. For ex- 
ample, a C2H2-vacancy for the structure in Fig. [I] (b) 
breaks the isotropy and can reduce C by ~ 12% in a spe- 
cific direction. Hydrogen frustratiorPED can a i so be a 
crucial type of defect, which would affect C, since the 
structure is locally compressed and Ao is influenced. 

We next consider the behavior of the system for higher 
values of the strain ranging from -0.02 to 0.45 in uniform 
expansion. For this purpose, we preferred a fully sym- 
metric hexagonal lattice with well defined high symmetry 
points in the BZ. Again the calculations are performed 
in a large (10x10) supercell as shown in Fig. [2] (a). The 
harmonic region can be taken between —0.02 < e < 0.02 
and it is followed by an anharmonic region where higher 
order terms are not negligible in the strain energy equa- 
tion. The anharmonic region is followed by a plastic re- 
gion where irreversible structural changes occur in the 
system and it transforms into a different structure after 
the yielding point. Fig. [2] (b) is the plot of strain en- 
ergy Es and its derivative (dE 5 (e) / de) with respect to 
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the applied strain. Two critical strain values can be de- 
duced from the plots. The first one, e ci , is the point 
where the derivative curve attains its maximum value 
and then starts to decrease. It occurs nearly at e = 0.23, 
where the C-C bond length is around 1.87 A. This means 
that for e > e Cl , the structure can be expanded under 
smaller tensions. The phonon frequencies, we calculated 
by using the force constant methocP^ are all positive 
throughout the BZ for e < e Cl , but the frequencies of 
longitudinal acoustic modes start to become imaginary 
for e > e ci , indicating an instability of 2D graphane un- 
der uniform expansion beyond e ei . Such phenomena is 
known as "phonon instability^ 18 * 19 !, where phonon fre- 
quencies fi n (k), get imaginary for specific wave vector k 
and branch index n. A detailed discussion can be found 
in Ref . |T5] and the references therein. Liu et a/P^ calcu- 
lated the critical strain values for graphene as 0.194 and 
0.266 for uniaxial tension in zigzag (x-) and armchair 
(y-) directions by using density functional perturbation 
theory (DFPT). 

The second critical point e C2 is the yielding point which 
is around e = 0.34. The C-C distance corresponding to 
e C2 is 2.02 A. Up to this point, the strain energy always 
increases and the system preserves its honeycomb-like 
structure. Upon the release of the tension, all the de- 
formation disappears and hence the system may return 
to its original size at e=0. Furthermore, the value of e C2 
is found to depend on various defects and the tempera- 
ture of the system. For H-vacancy, we found that e C2 is 
lowered to ~0.21 as shown in Fig. [2] (c). As for C+H- 
vacancy, which corresponds to a hole at one corner of 
hexagon, e C2 is further lowered to 0.13 as shown in Fig. [2] 
(d). We also examined the effect of ambient tempera- 
ture on the yielding strain. Ab-initio molecular dynamic 
calculations (lasting 2 ps with time steps of 2xl0~ 15 sec- 
onds) indicate that e C2 =0.34 corresponding to T=0 K 
is reduced to 0.20 at T=300 K and is further reduced 
to 0.18 at T=600 K. Apparently, the yielding of perfect 
graphane under uniform strain at e C2 can only occur for 
ideal conditions. For e Cl < e < e C2 the system is in a 
meta-stable state. The long wave length perturbations, 
vacancy defects, as well as high temperature effects lead 
e C2 decrease to the strain values around e Cl . After the 
yielding point, where e > e C2 , the plastic range sets in 
with irreversible deformations. This range, however, is 
beyond the scope of this paper. 

We finally investigate the variation of the electronic 
properties of graphane with the uniform strain. The ef- 
fect of strain on the buckling is found to be minute. It 
decreases from 0.46 A to 0.43 A as e increases from 
to 0.30. Also, C-H bonds are shortened only 1% in this 
range of strain. The binding energy of a single hydro- 
gen in (10x10) supercell increases from 4.79 to 5.02 up 
to e=0.20. Normally, graphane is a semiconductor with 
a wide direct band gap of 3.54 eV calculated by DFT- 
GGA, but our calculations show that this gap can in- 
crease to 5.66 eV after G W corrections. On the other 
hand, recent GW (5.97 eV^ and GW (5.4 eV)P cor- 
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Figure 3: (Color Online) The variation of energy band gaps 
with (2D) uniform expansion. The band gaps obtained both 
from GGA (green triangles) and G W (blue squares) calcu- 
lations increase with increasing strain up to e = 0.15, passes 
through a maximum, then decrease until the yielding point. 
The band gaps are given on the left and the strain energies 
are given on the right. Three panels show how the bands at 
the edge of conduction band are modified with strain. 



rections report slightly different values depending on the 
method and parameters used. More recently DFT-LDA 
calculations^ found the band gap as 3.6 eV. Figure. [3] 
shows the variation of GGA and Go Wo band gap val- 
ues with respect to the strain for uniform expansion in 
the elastic region. While the lowest conduction band is 
raised with strain in the first and second panels; in the 
third panel, the second ir* band is lowered steadily and 
dips in the gap for e > 0.15. Dramatic variation of the 
band gap with the strain suggests that graphane can be 
used as a strain gauge at nanoscale. 

In summary, we revealed the elastic constants of 
graphane indicating that it has a quite high in-plane 
stiffness and very low, perhaps the lowest Poisson's ra- 
tio among known monolayer honeycomb structures. We 
showed that the band gap of graphane can be modified 
significantly by applied strain in the elastic range. It is 
suggested that elastic deformation can be used for further 
functionalization of graphane and hence for monitoring 
its chemical and electronic properties. 

Part of the computations have been provided by 
UYBHM at Istanbul Technical University through a 
Grant No. 2-024-2007. 
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